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The ultra-fast magnetisation relaxation rates during the laser-induced magnetisation process are analyzed 
in terms of the Landau-Lifshitz-Bloch (LLB) equation for different values of spin S. The LLB equation is 
equivalent in the limit S — > °° to the atomistic Landau-Lifshitz-Gilbert (LLG) Langevin dynamics and for 5=1/2 
to the M3TM model [B. Koopmans, et al. Nature Mat. 9 (2010) 259]. Within the LLB model the ultra-fast 
demagnetisation time (xm) and the transverse damping (oc^) are parameterized by the intrinsic coupling-to-the- 
bath parameter X, defined by microscopic spin-flip rate. We show that for the phonon-mediated Elliott- Yafet 
mechanism, X is proportional to the ratio between the non-equilibrium phonon and electron temperatures. We 
investigate the influence of the finite spin number and the scattering rate parameter X on the magnetisation 
relaxation rates. The relation between the fs demagnetisation rate and the LLG damping, provided by the LLB 
theory, is checked basing on the available experimental data. A good agreement is obtained for Ni, Co and 
Gd favoring the idea that the same intrinsic scattering process is acting on the femtosecond and nanosecond 
timescale. 
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I. INTRODUCTION 

Magnetisation precession and the spin-phonon relaxation 
rates at picosecond timescale were considered to be the limit- 
ing factor for the speed of the magnetisation switching^, un- 
til using optical excitation with fs pulsed lasers the possibil- 
ity to influence the magnetisation on femtosecond timescale 
was demonstrated 3 " 6 . The ultra-fast laser-induced demagneti- 
sation immediately became a hot topic of solid state physics 
due to an appealing possibility to push further the limits of 
operation of magnetic devices. This ultra-fast process has 
now been shown to proceed with several important charac- 
teristic timescales 6 : (i) the femtosecond demagnetisation with 
timescale %m (ii) the picosecond recovery with timescale Xe 
and (iii) the hundred picoseconds -nanosecond magnetisation 
precession, traditionally characterized by the ferromagnetic 
resonance frequency (£>fmr and the Landau-Lifshitz-Gilbert 
damping parameter (Xllg ( see Fig.l). 

The physics of the magnetisation changes on femto-second 
timescales is obviously not-trivial and will require novel the- 
ories within the relativistic quantum electrodynamics of many 
electron systems. From theoretical point of view, the existing 
models try to answer an open question of the role of differ- 
ent subsystems (photons, phonons, electrons and spins) in the 
ultra-fast angular momentum transfer^. This common goal is 
stimulated by experimental findings provided by the XMCD 
measurements showing the important role of the spin-orbit 
interactions 8 -* 9 -. For the present state of art quantum mechani- 
cal descriptions^^"— of ultra-fast demagnetisation processes 
involve unavoidable simplifications and sometimes even some 
ad-hoc assumptions necessary to explain experimental find- 
ings, such as reduced exchange interactions, enhanced spin- 
orbit coupling or a Gaussian distribution of occupied states 
around the Fermi level. While some degree of agreement has 
been achieved in modelling of the ultra-fast demagnetisation 
(Xm) scaled, the modelling of all three ultra-fast demagnetisa- 
tion rates within the same approach is outside the possibilities 
of the quantum mechanical approaches. 

The three-temperature (3T) phenomenological model in- 



volves the rate equations for the electron, phonon and spin 
temperatures (energies) 10 ' 15 " —. Recently, it has been shown 
that the introduction of the spin temperature is not adequate-^ 
since the spin system is not in the equilibrium on the fem- 
tosecond timescale. It has been suggested to couple the spin 
dynamics to the two-temperature (2T) model for phonon and 
electron temperatures^—. These models are based on the en- 
ergy flow picture and leave unidentified the angular momen- 
tum transfer mechanism and the underlying quantum mech- 
anism responsible for the spin flip^ 2 .. They essentially inter- 
pret the ultra-fast demagnetisation as "thermal" processes, un- 
derstanding the temperature as energy input from photon to 
electron and then to the spin system. By using these mod- 
els the important role of the linear reversal path in the femto- 
second demagnetisation has been identifie d 23 ' 24 . The compar- 
ison with experiment seems to indicate that in order to have 
magnetisation switching in the ultra-fast timescale, a com- 
bined action of "heat" and large field coming from the inverse 
Faraday effect is necessary^ 4 -. 

The most successful recent phenomenological models de- 
scribing the ultra-fast magnetisation dynamics are (i) the 
Langevin dynamics based on the Landau-Lifshitz-Gilbert 
(LLG) equation and classical Heisenberg Hamiltonian for lo- 
calized atomic spin moment s 18 ' 19 , (ii) the Landau-Lifshitz- 
Bloch (LLB) micromagnetic a 21 ' 22 and (iii) the Koopmans's 
magnetisation dynamics model (M3TM) 25 - . The spin dy- 
namics could be coupled to the electron temperature from the 
2T model, underlying the electronic origin of the spin-flip 
proces s 18 ' 19 ' 21 ' 22 ' 24 or to both electron and phonon temper- 
atures, underlying the Elliott- Yafet mechanism mediated by 
phonons 25 . When the 2T model was carefully parameterised 
from the measured reflectivity, it gave an excellent agreement 
with the experiment in Ni 22 using the former approach or in 
Ni, Co and Gd using the latter approach^ 5 -. 

In the classical derivation of the LLB equation the ther- 
mal averaging has been performed analytically within the 
mean field (MFA) approximation 26 . Thus, the LLB equa- 
tion for classical spins (S — > °°) is equivalent to an ensem- 
ble of exchange-coupled atomistic spins modelled by stochas- 



tic LLG equation a 20 ' 27 . At the same time, in some cases the 
LLB equation may be preferable with respect to the atomistic 
Heisenberg model, since being micromagnetic it can incorpo- 
rate quantum nature of magnetism and the quantum deriva- 
tion of LLB also exists 2 ^. In particular the limits of validity 
for the statistical mechanics based on the classical Heisenberg 
model for the description of materials with delocalized mag- 
netism of ^-electrons in transition metals or magnetism of/- 
electrons in rare earths are not clear. An alternative statistical 
simplified description of c/-metals consists of a two level sys- 
tem with spin-up and spin-down bands (i.e. S = ±1/2), as 
has been done by B. Koopmans et al^. Their model, as we 
show in the present article, is also equivalent to the quantum 
LLB equation with spin S = 1/2. An additional advantage 
in the use of the LLB equation is the possibility to model 
larger spatial scale s - 20 ' 21 . Therefore the LLB micromagnet- 
ics is an important paradigm within the multiscale magnetisa- 
tion dynamics description. The LLB equation has been shown 
to describe correctly the three stages of the ultra-fast demag- 
netisation processes: the sub-picosecond demagnetisation, the 
picosecond magnetisation recovery and the nanosecond mag- 
netisation precession 2 ^"—, see FigQ] 

The intrinsic quantum mechanical mechanisms responsi- 
ble for the ultra-fast demagnetisation in the LLB model 
are included in the intrinsic coupling-to-the-bath parame- 
ter Ai 22 i 28 . The coupling process is defined by the rate of 
the spin flip. Several possible underlying quantum mech- 
anisms are currently under debate: the Elliott- Yafet (EY) 
electron scattering mediated by phonons or impuritie s - 13 ' 25 , 
or other electrons^ 4 , and electron-electron inelastic exchange 
scatterin g 29 ' 30 . By combining the macroscopic demagnetisa- 
tion equation (M3TM model) with the rate of spin flip calcu- 
lated on the basis of full Hamiltonian, Koopmans et alJ£- have 
been able to relate the ultra-fast demagnetisation time Zm with 
the spin flip rate of the phonon-mediated Elliott- Yafet scatter- 
ing. The authors fitted experimental demagnetisation rates in 
Ni, Co, Gd to the phenomenological M3TM model and found 
them to be consistent with the values estimated on the basis of 
ab-initio theory^. The coupling-to-the-bath parameter A, (mi- 
croscopic damping parameter in atomistic LLG model) should 
be distinguished from that of the macroscopic damping O-llg 
(0Cx in the LLB model), a more complicated quantity which 
includes the magnon-magnon processes. 

The first attempt to relate the sub-picosecond demagnetisa- 
tion time with the macroscopic damping processes was given 
by Koopmans et alk who suggested the relation 1m ~ 1 / Q-llg- 
Subsequently and with the aim to check this relation several 
experiments in doped permalloy were performed 3 - 2 - - — . The 
permalloy thin films were doped with rare earth impurities, al- 
lowing to increase in a controlled way the damping parameter 
OCllg- The effect on the demagnetisation time %m was shown 
to be opposite 3 - 4 - or null 3 - 2 -, in contrast to the above relation. 
However, it should be noted that the analysis leading to this 
expression was performed in terms of the Landau-Lifshitz- 
Gilbert equation, relating the ultra-fast demagnetisation time 
1m to the transverse damping without taking into account their 
temperature dependence. Moreover, one should take into ac- 
count that the rare-earth impurities may introduce a different 
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Figure 1. Characteristic time scales in ultrafast laser-induced mag- 
netisation dynamics experiments. The curve is obtained by the in- 
tegration of the Landau-Lifshitz-Bloch equation coupled to the two- 
temperature model with the parameters from Ref|2ll For the mod- 
elling of precession the applied field H ap = IT at 30 degrees was 
used. 



scattering mechanism with a slower timescale 3 - 3 -. 

Partially basing on the above mentioned experimental re- 
sults and from a general point of view, the longitudinal relax- 
ation (the ultra-fast demagnetisation rate Tm) and the trans- 
verse relaxation (the LLG damping olllg) ma y be thought 
to be independent quantities. Indeed, different intrinsic and 
extrinsic mechanisms can contribute to the demagnetisation 
rates at different timescales. One can, for example, men- 
tion that during the femtosecond demagnetisation the electron 
temperature is often raised up to the Curie temperatur e) 22 ' 24 . 
At this moment, the high frequency THz spinwave o 35 ' 36 in- 
cluding the Stoner excitations^ contribute. At the same time, 
the transverse relaxation is related to the homogeneous pre- 
cessional mode. The LLB equation takes care of the different 
natures of longitudinal and transverse relaxation, arising from 
the spin disordering. The LLB model calculates them inde- 
pendently but basing on the same intrinsic scattering mecha- 
nism parameterized by the parameter X. The increment of the 
number of scattering events is mimicked by the increases of 
the electron temperature. Consequently, the relation between 
the ultra-fast demagnetisation and precession remains valid 
but with a temperature-dependent correction. If this relation is 
confirmed experimentally, a unique intrinsic coupling param- 
eter means that the same main microscopic mechanism is act- 
ing on both timescales. In the present article we will show that 
the analysis of the available experimental data seems to indi- 
cate towards this possibility, at least in pure transition metals 
such as Ni or Co and in rare earth metal Gd. We did not find 
validity of the corresponding relation in Fe. 

Up to now only classical version (S — > °°) of the LLB 
equation was used to model the ultra-fast demagnetisation 
processes 



In the present article we show the impor- 



tant role of the choice of the quantum spin value, resulting 
in the differences in the corresponding longitudinal relaxation 
times. The article is organized as follows. In section II we 
present different formulations of the quantum LLB model and 
its main features for different spin values S. In section III 
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we present results on the modelling of the demagnetisation 
processes within LLB models with different choices of the 
quantum spins number S and of the intrinsic scattering mech- 
anisms. In section IV we present our attempts to link the 
ultra-fast demagnetisation rates in transition metals and Gd 
and comparison with available experimental data. Section 
V concludes the article. In the Appendix to the article we 
demonstrate the equivalence of the LLB model with 5=1/2 
and the M3TM model by B. Koopmans et alJ^. 

II. THE LANDAU-LIFSHITZ-BLOCH MODEL WITH 
QUANTUM SPIN NUMBER 5. 

The LLB equation for a quantum spin was derived from the 
density matrix approach^. Although the model Hamiltonian 
was rather the simplest form of the spin-phonon interaction, 
the generalization of the approach should be possible to more 
complex situations. The macroscopic equation for the mag- 
netisation dynamics, valid at all temperatures, is written in the 
following form: 



with qs = 3Tcm e / [2(S + l)T], where S is the quantum spin 
number and Tc is the Curie temperature. In the case S —> °° 
the damping coefficients have the forms used in several previ- 
ously published works- , suitable for the comparison with the 
Langevin dynamics simulations based on the classical Heisen- 
berg Hamiltonian and in agreement with the m 20 ! 27 . 

Eq.(Q~|i is singular for T > Tc, in this case it is more con- 
venient to use the LLB equation in terms of the variable 
m = M/M e (T = OK) 21 . The corresponding LLB equation is 
indistinguishable from Eq.(fl]i but with different relaxational 
parameters fi|| = m e ap fix = m e a^ and fix = fi|| for T > Tc, 
in this case the contribution of temperature to H e g- [the sec- 
ond term in Eq.©] is (- 1 /%\\ )[1 - 3T c m 2 /5(T - T c )m]m. Al- 
though this formulation is more suitable for the modelling of 
the laser-induced demagnetisation process, during which the 
electronic temperature is usually raised higher than Tc, it is 
the expression (O which should be compared with the trans- 
verse relaxation parameter (Xllg due to the similarity of the 
formulation of the Eq.(Q]) with the macromagnetic LLG equa- 
tion. In the classical case and far from the Curie temperature 
T <c T c , X = ocx = fi ± (a LLG ). 



n = y[nxH] + — ^[n-H e ff]n r nx [nxH eff 
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where n = M/M e (T) = m/m e is the reduced magnetisation, 
normalized to the equilibrium value M e at given temperature T 
and m = M/M e (T = OK). The effective field H e ff, contains all 
usual micromagnetic contributions, denoted by Hj nt (Zeeman, 
anisotropy, exchange and magnetostatic) and is augmented by 
the contribution coming from the temperature 
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where %\\{T) = (dm/dH)H^o is the longitudinal susceptibil- 
ity . The LLB equation contains two relaxational parame- 
ters: transverse ax and longitudinal OCm, related to the intrinsic 
coupling-to-the-bath parameter X. In the quantum description 
the coupling parameter X contains the matrix elements repre- 
senting the scattering events and, thus, is proportional to the 
spin-flip rate due to the interaction with the environment. This 
parameter, in turn, could be temperature dependent and, in our 
opinion, it is this microscopic parameter which should be re- 
lated to the Gilbert parameter calculated through ab-initio cal- 
culations as in Refsi 38,39 , since the contribution coming from 
the spin disordering is not properly taken into account in these 
models. In the quantum case the temperature dependence of 
the LLB damping parameters is given by the following ex- 
pressions: 
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Figure 2. (Up) The transverse damping parameter ax (Cllg) as 
a function of temperature within the LLB model for different spin 
values S. The intrinsic coupling parameter was set to X = 0.03. 
(Down) The longitudinal relaxation time Tn as a function of tem- 
perature within the LLB model for different spin values S. The 
temperature-dependent magnetisation and the longitudinal suscepti- 
bility X|| were evaluated in both cases in the MFA approach using the 
Brillouin function. 



In the "thermal" model the nature of the longitudinal and 
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the transverse relaxation differs from the point of view of 
characteristic spinwave frequencies. The transverse relaxation 
(known as the LLG damping) is basically the relaxation of the 
FMR mode. The contribution of other spinwave modes is re- 
duced to the thermal averaging of the micromagnetic parame- 
ters and the main effect comes from the decrease of the mag- 
netisation at high temperature. Consequently, the transverse 
damping parameter increases with temperature (see Figf2]i, 
consistent with atomistic modelling results 2 ^, and well-known 
FMR experiment a 37 ' 41 . 

On the contrary, the main contribution to the longitudinal 
relaxation comes from the high-frequency spin waves. This 
process occurs in a strong exchange field. As a result, the lon- 
gitudinal relaxation time (the inverse longitudinal relaxation) 
is much faster and increases with temperature, known as crit- 
ical slowing down, see Figj2] This slowing down has been 
shown to be responsible for the slowing down of the femto- 
second demagnetisation time 1m as a function of laser pump 
fluenc y 18 ' 22 . The characteristic longitudinal timescale is not 
only defined by the longitudinal damping parameter ([3]l but 
also by the temperature-dependent longitudinal susceptibility 
%\\ (7")—, according to the following equation: 



xii (n 
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As it can be observed in Fig. [2] the transverse relaxation 
parameter OCx(oCzxg) an d the longitudinal relaxation time Tii 
have a strong dependence on the quantum spin number S cho- 
sen to describe system's statistics. We conclude here about 
the occurrence of quite different relaxation rates for the two 
extreme cases 5=1/2 and S = °°. 

B. Koopmans et al. recently used a different equation 
to describe the ultrafast demagnetisation dynamics 2 ^,, called 
M3TM model: 



dm T„ ( ( mTc 

= Rm— 1 — mcoth 

dt T c \ \ T e 



(6) 



Eq.© has been obtained through the general Master equation 
approach for the dynamics of the populations of a two level 
system (spin 5=1/2 was used) with the switching probabil- 
ity evaluated quantum-mechanically for the phonon-mediated 
EY spin-flips. Here T p and T e are phonon and electron tem- 
peratures, respectively, and R is a material specific parameter, 
related to the spin-flip probability in the phonon-mediated EY 
scattering events a s f, as 
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where V a and /j at are the atomic volume and magnetic moment, 
respectively, G ep is the electron-phonon coupling constant, kg 
is the Boltzmann constant, /jg is the Bohr magneton and Ed is 
the Debye energy. This equation has allowed to fit the ultra- 
fast demagnetisation time (Tm) obtaining the values of R in 
Ni, Co and Gd 2 ^ and relating them to the phonon-mediated 
EY scattering rates a s f. 

As we show in the Appendix, the M3TM equation (O cor- 
responds to the longitudinal part of the LLB equation with 



thermal field only (H; nt = 0) and with spin S = 1/2, i.e. it is 
equivalent to 



dm ^ 

— =ya\\H eS . (8) 
dt " 

This gives a relation between the intrinsic coupling parameter 
X and the material specific parameter R and finally with the 
phonon-mediated EY spin-flip probability a s f via the formula: 



3R /J a t T p _ T p 



(9) 



Thus the two approaches are reconciled, provided that the 
temperature-dependent coupling rate (O is used in the LLB 
equation, in contrast to other work a ' ■ where the coupling 
X is considered to be temperature-independent. Combining 
expressions <(5j (|7]i and (O, one can immediately see that in 
the case of the phonon-mediated EY process, the longitudinal 
relaxation time is determined by 



X|| 
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In Ref and basing on the phonon-mediated EY picture, the 
classification of materials on the basis of the "magnetic in- 
teraction strength" parameter fi a t/J was proposed, where J is 
the material exchange parameter. According to the expression 
above, the demagnetisation rate depends on more parameters, 
among which the important one is also the electron-phonon 
coupling G ep defining how fast the electron system can pass 
the energy to the phonon one. Another important parameter is 
the microscopic spin-flip rate a s f. Comparing to the B. Koop- 
mans et alJ£- materials classification, the longitudinal suscep- 
tibility in Eq.JTOt is indeed defined by the value of the atomic 
moment /j at and by the fact that this function rapidly increases 
with temperature and diverges close to Tq « /. At T ~ Tc one 
obtains a simple linear relation 27 X|| x aW-A thus showing the 
dependence of the demagnetisation rate on this parameter, as 
suggested in Refl25i 

In the case of the phonon-mediated EY process the tem- 
perature dependence of the longitudinal relaxation is coming 
from the longitudinal susceptibiliy only (cf. Eq. (fTUli). as op- 
posed to the case X = const (cf. Eq.©). (We do not discuss 
here the possibility that the phonon-mediated EY spin-flip rate 
a s f may be also temperature dependent.) However, the tem- 
perature dependence of the susceptibility is characterized by 
its exponential divergence close to Tc. In these circumstances 
an additional linear temperature dependence provided by the 
longitudinal damping is difficult to distinguish in the fitting 
procedure of experimental data. 



in. MODELLING OF THE LASER-INDUCED 
ULTRA-FAST DEMAGNETISATION WITHIN THE LLB 
MODELS. 

In the spirit of Refs i 18 i 20 " 22 ' 25 for the modelling of ultra-fast 
demagnetisation dynamics, the LLB equation may be coupled 
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to the electron temperature T e only, understanding the elec- 
trons as the main source for the spin-flip mechanism^™* — 
or to both phonon and electron temperatures in the spirit of 
the phonon-mediated Elliott- Yafet process 2 - 5 -. In both cases it 
is the electron temperature T = T e which couples to the mag- 
netisation in the LLB formalism, since the phonon tempera- 
ture could only enter into the temperature dependence of the 
coupling-to-the bath parameter X via Eq.© . Note that the 
temperature T is not the spin temperature, since the resulting 
dynamics is taking place out-of-equilibrium. 

The electron T e and phonon T p temperatures are taken from 
the two-temperature (2T) mode l 15 ' 45 ' 46 . Within this model 
their dynamics is described by two differential equations: 

C e -^- = -G ep (T e -T p )+P(t), 

C p ^=G ep (T e -T p ). (11) 

Here C e = y e T e (y e = const) and C p are the specific heats of 
the electrons and the lattice. The Gaussian source term P(t) is 
a function which describes the laser power density absorbed 
in the material. The function P(t) is assumed to be propor- 
tional to the laser fluence F with the proportionality coeffi- 
cient which could be obtained from the long time scale de- 
magnetization data (for which T e = T p )*£. The dynamics of 
the electron temperature can be also measured directly in the 
time-resolved photoemission experiment 4 ^. 

The first of Eqs. (fTTT > may also include a diffusion term 
V-(KV z r c ) taking into account a final penetration depth of 
the deposited energy into the film thickness^ 5 - and a term, 
C e (T e — 3QQK)/z t i 1 describing the heat diffusion to the exter- 
nal spaced. In the present article, the parameters for the 2T- 
model were taken either from Koopmans et alM- or from U. 
Atxitia et al.~ (for Ni only), where they were carefully pa- 
rameterized through the reflectivity measurements. The Ni 
(Co, Gd etc) parameters, such as magnetisation as a function 
of temperature were taken assuming the Brilloiun (Langevin 
for S — > oo) function. 

The coupling of the 2T model to the LLB equation ade- 
quately describes all three stages of the ultra-fast demagneti- 
sation rates: sub-ps demagnetisation, ps recovery and sub-ns 
precessio n 21 ' 22 , see Fig. 1 . As a consequence of the temper- 
ature dependence of both longitudinal damping and suscepti- 
bility, and since the temperature is dynamically changed ac- 
cording to Eqs.dTTt. the longitudinal relaxation time is time- 
dependent via Eq.©. It is also strongly dependent on the 
parameters of the 2T model and its dynamics is not simple. 
Consequently, the sub-ps ultra-fast demagnetisation generally 
speaking is not exponential and cannot be described in terms 
of one relaxation time Xm- Simple analytical expression is 
possible to obtain with the supposition of a square-shaped 
temperature pulse 23 . The two-exponential fitting is also often 
use d 22 ' 36 . In our approach the fs demagnetisation is fitted di- 
rectly to the solution of the LLB equation without assumption 
of the one- or two-exponential decay. However, to comply 
with the existing approaches, we still discuss the demagneti- 
sation rate in terms of a unique parameter Xm- 

In the experiment performed in the same material the only 



remaining fitting parameter for the LLB model is the coupling 
parameter X. The choice of X together with the parameters 
of the 2T model defines all magnetisation rates. In Figf3]we 
present modelling of the ultra-fast demagnetisation and re- 
magnetisation for various values of the coupling parameter X, 
chosen to be independent on temperature, as in Ref. [22[ If 
for some reason the scattering channel was suppressed, this 
would lead to a small scattering rate and consequently a small 
demagnetisation and a slow recovery. Indeed, the value of X 
for Gd was found to be 60 times smaller than for Ni (see Table 
I). This small value of X assures a large delay in the magnetis- 
arion relaxation towards the equilibrium electron temperature. 
Thus this parameter defines the diversity of the demagnetisa- 
tion rates in larger extend than the ratio /J at /J, suggested in 
RefEBI and discussed in the previous subsection. 
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Figure 3. The result of integration of the LLB model (5 — > °°) with 
different parameters X (increasing from top to the bottom). In this 
case the the 2T model parameters were taken from Refill with laser 
fluence F = 30 mj/cm 2 

Another parameter strongly influencing the demagnetisa- 
tion rates is the phonon-electron coupling G ep defining the 
rate of the electron temperature equilibration time. This is 
the main parameter governing the magnetisation recovering 
time Te ■ Indeed, in Ref.— the phonon-electron coupling G ep 
was chosen to be 20 times smaller for Gd than for Ni. By 
adjusting this parameter, the ultra-slow demagnetisation rates 
observed in TbFe alloy 4 -^, GdM and in half-metals 5 ^ as well 
as the two time-scales demagnetisatio n 25 ' 49 are also well- 
reproduced (see, as an example, FiglUi. Within this model 
the two-time scale process consists of a relatively fast demag- 
netisation (however much slower than in Ni), defined by the 
electron temperature and small value of X, followed by a much 
slower process due to a slow energy transfer from the electron 
to the lattice system. 

As it was mentioned in the previous subsection, the 
phonon-mediated EY mechanism predicts the coupling to the 
bath parameter X to be dependent on the ratio between the 
phonon and electron temperature through the relation (O. A 
decrease of X up to two times at high fluencies is observed 
for Ni and Co. The analysis of the data presented in Ref. [25] 
and|47|for Gd has shown that during the demagnetisation pro- 
cess the ratio T e /T p has increased almost 6 times. In Fig(5] 
we present the magnetisation dynamics for Ni evaluated for 
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Figure 4. The result of integration of the LLB model (5 — > °°) with 
constant Xq = 0.0015 (see Table I). In this case the 2T model param- 
eters were taken from Ref|25l corresponding to Gd. 



two laser pulse fluencies, assuming various values of the spin 
S and temperature-dependent and independent A, values. Note 
quite different demagnetisation rates at high fluency for two 
limiting cases S — 1/2, used in Refl25l and S — °°, used in 
Refl22l The differences in the choice of X are pronounced 
at high pump fluency but are not seen at low fluency. One 
can also hope that in the fitting procedure of experimental 
data it would be possible to distinguish the two situations. 
Unfortunately, the fitting to experimental data procedure is 
complicated and the changes coming from the two cases de- 
scribed above are competing with several different possibil- 
ities such as an additional temperature dependency in elec- 
tron or phonon specific heats^i. Additionally, we would like 
to mention different electron-phonon coupling constants G ep 
used in Refs. |22| and Fitting to experimental data from 
Refl25l for Ni for high fluence, we have found that the case 
of the temperature-dependent X = Xo(T p /T e ) can be equally 
fitted with the temperature-independent A, « Xq/2. To answer 
definitely which fitting is better, more experimental data pro- 
moting one or another intrinsic mechanism and varying laser 
fluency is necessary. 



IV. LINKING DIFFERENT TIMESCALES 

Since the longitudinal relaxation occurs under strong ex- 
change field and the transverse relaxation - under external 
applied field, their characteristic timescales are quite differ- 
ent. However, the LLB equation provides a relation be- 
tween the ultra-fast demagnetisation (longitudinal relaxation) 
and the transverse relaxation (ordinary LLG damping param- 
eter) via the parameter A,o (X = Xq or X = Xo(T p /T e ) for 
T p = T e ). The two demagnetisation rates could be measured 
independently by means of the ultra-fast laser pump-probe 
technique^. It has been recently demonstrated 5 - that the 
damping of the laser-induced precession coincides with the 
measured by FMR in transition metals. By separate measure- 
ments of the two magnetisation rates, the relations (|4|i and (0 
given by the LLB theory could be checked. This can pro- 
vide the validation of the LLB model, as well as the answer to 
the question if the same microscopic mechanism is acting on 




-0.5 



t [ps] 

Figure 5. Magnetisation dynamics during laser-induced demagneti- 
sation process calculated within the LLB model with different spin 
numbers and for two laser- fluencies F = 10 mJ/cm 2 (upper curves) 
and F = 40 mJ/cm 2 (bottom curves). Ni parameters from Refill 
were used. The symbols are calculated with the LLB equation with 
the intrinsic damping parameter using a constant Xq = 0.003 value, 
and the solid lines with the LLB equation and the intrinsic coupling 
with the temperature dependent X = Xq (Tp/T e ) . 



femtosecond and picosecond timescales. Unfortunately, the 
damping problem in ferromagnetic materials is very compli- 
cated and the literature reveals the diversity of measured val- 
ues in the same material, depending on the preparation condi- 
tions. 

Thus, to have a definite answer the measurement on 
the same sample is highly desired. The measurements of 
both ax and 1m are available for Ni^ 2 . where an excel- 
lent agreement between ultra-fast magnetisation rates via a 
unique temperature-independent parameter X = 0.04 has been 
reported 2 ^. The results of the systematic measurements of 1m 
are also available for Ni, Co, Gd in Ref. |25[ as well as for 
Fa^. The next problem which we encounter here is that the 
demagnetisation rates strongly depend on the spin value S, as 
is indicated in Figs. |2]and|5] The fitting of experimental data 
using LLB model with different S values results in different 
values of the coupling parameter Xq. The use of S = 1/2 
valued or S = °° valued is quite arbitrary and these values 
do not coincide with the atomic spin numbers of Ni,Co, Gd. 
Generally speaking, for metals the spin value is not a good 
quantum number. The measured temperature dependence of 
magnetisation, however, is well fitted by the Brillouin func- 
tion with 5=1/2 for Ni and Co and S = 7 /2 for Gd^i These 
are the values of S which we use in Table I. 

Consequently in Table I we present data for the coupling 
parameter Xq extracted from Ref 25. Differently to this article, 
for Gd we corrected the value of the parameter R to account 
for a different spin value by the ratio of the factors, i.e. R Sl = 
(/^//sj/?* 2 with 



fs 



2qs 



l 



smh(2q s ) m 2 s xf ' 



(12) 



where the parameters are evaluated at 120K using the MFA 
expressions for each spin value S. The data are evaluated 
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Material S RsZ Xq rx ± (Xllg 

Ni 1/2 17.2 0.0974 0.032 0.019 42 -0.0282i 
Co 1/2 25.3 0.179 0.025 0.0036^-0.00622-0.01 1-£ 

Gd 7/2 0.009 0.0015 0.00036 0.0005*2 

Table I. The data for ultra-fast demagnetisation rate parameters for 
three different metals from ultrafast demagnetization rates and from 
FMR mesurements. The third column presents the demagnetisation 
parameter R from Ref. [25L corrected in the case of Gd for spin 
S = 7/2. The fourth column presents the value of the Xq parame- 
ter, as estimated from the M3TM models^ and the formula Eq.(|9). 
The fifth column presents the data for ax estimated via the LLB 
model Eq.«(4j and the Xq value from the third column, at room tem- 
perature T = 300A" for Co and Ni and at T = 120A" for Gd . The 
last column presents the experimentally measured Gilbert damping 
collected from different references. 

for the phonon-mediated EY process with the temperature- 
dependent parameter X via the expression (0. The value of 
the Gilbert damping parameter oc^ was then estimated through 
formula © at 300K (for Ni and Co) and at 120K for Gd. Note 
that for temperature-independent X = Xq the resulting A-o and 
OCx values are approximately two times smaller for Ni and Co. 
The last column presents experimental values for the same pa- 
rameter found in literature for comparison with the ones in the 
fifth column, estimated through measurements of the ultra- 
fast demagnetisation times Xm and the relation provided by 
the LLB equation. 

Given the complexity of the problem, the results presented 
in Table I demonstrate quite a satisfactory agreement between 
the values, extracted from the ultra-fast demagnetisation time 
Zm and the Gilbert damping parameter via one unique 
coupling-to-the-bath parameter X. The agreement is particu- 
larly good for Ni, indicating that the same spin flip mechanism 
is acting on both timescales. This is true for both experiments 
in Refsl22l and HH For Co the value is some larger. For the 
temperature-independent X, the resulting value is two times 
smaller and the agreement is again satisfactory. We would 
like to note that no good agreement was obtained for Fe. The 
reported damping values^ are 5-10 times smaller as estimated 
from the demagnetisation rates measured in Ref. [HI 

V. CONCLUSIONS 

The Landau-Lifshitz-Bloch (LLB) equation provides a mi- 
cromagnetic tool for the phenomenological modelling of the 
ultra-fast demagnetisation processes. Within this model one 
can describe the temperature-dependent magnetisation dy- 
namics at arbitrary temperature, including close and above 
the Curie temperature. The micromagnetic formulation can 
take into account the quantum spin number. The LLB model 
includes the dynamics governed by both the atomistic LLG 
model and the M3TM model by Koopmans et al^. In 
the future it represents a real possibility for the multiscale 
modelling^. 

We have shown that within this model the ultra-fast de- 
magnetisation rates could be parameterized through a unique 



temperature dependent or independent parameter X, defined 
by the intrinsic spin-flip rate. The magnetisation dynamics 
is coupled to the electron temperature through this parameter 
and is always delayed in time. The observed delay is higher 
for higher electron temperature. This is in agreement with the 
experimental observation that different materials demagnetize 
at different rate a 25 i 50 and that the process is slowed down with 
the increase of laser fluency. We have shown that for the 
phonon-mediated EY mechanism the intrinsic parameter A, is 
dependent on the ratio between phonon and electron temper- 
atures and therefore is temperature dependent on the femto 
second - several picosecond timescale. The LLB equation can 
reproduce slow demagnetizing rates observed in several ma- 
terials such as Gd, TbFe and half metals. This is in agree- 
ment with both phonon-mediated EY picture since in Gd a 
lower spin-flip probability was predicted and also with the in- 
elastic electron scattering picture, since the electron diffusive 
processes are suppressed in insulators and half-metal a 31 ' 50 . 
However, we also stress the importance of other parameters 
determining the ultra-fast demagnetisation rates, such as the 
electron-lattice coupling. 

The macroscopic damping parameters (longitudinal and 
transverse) have different natures in terms of the involved 
spinwaves and in terms of the timescales. Their temperature 
dependence is different, however, they are related by the spin- 
flip rate. We have tried to check this relation in several tran- 
sition metals such as Ni, Co, Fe and the rare-earth metal Gd. 
A good agreement is obtained in Co and Gd and an excel- 
lent agreement in Ni. This indicates that on both timescales 
the same main microscopic mechanism is acting. In Ni the 
agreement is good both within the assumptions X = Xo and 
X = XoT p /T e . In Co the agreement seems to be better with the 
temperature-independent parameter X = Xo which does not in- 
dicate towards the phonon-mediated EY mechanism. How- 
ever, given a small discrepancy and the complexity of the 
damping problem, this conclusion cannot be considered defi- 
nite. We can neither exclude an additional temperature depen- 
dence of the intrinsic scattering probability (i.e. the parame- 
ter A-o) for both phonon-mediated EY and exchange scattering 
mechanisms which was not taken into account. 

An open question is the problem of doped permalloy where 
an attempt to systematically change the damping parameter by 
doping with rare-earth impurities was undertaken 3 - 3 - in order to 
clarify the relation between the LLG damping and the ultra- 
fast demagnetisation rat o 32 i 34 . The results are not in agree- 
ment with the LLB model. However in this case we think 
that the hypothesis of the slow relaxing impurities presented 
in Refl34l might be a plausible explanation. Indeed, if the 
relaxation time of the rare earth impurities is high, the stan- 
dard LLB model is not valid since it assumes an uncorrected 
thermal bath. The correlation time could be introduced in 
the classical spin dynamics via the Landau-Lifshitz-Miyasaki- 
Seki approach^. It has been shown that the correlation time 
of the order of 10 fs slows down the longitudinal relaxation 
independently on the transverse relaxation. Thus in this case, 
the modification of the original LLB model to account for the 
colored noise is necessary. 
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We expand around equilibrium m e = 1/2 (<?<?) the small 
quantity 1 — Bx^/m 



(A4) 



Appendix A 



where 8m = m — m e . Next, we expand m around to 2 



To show the equivalence between the LLB model with 
5=1/2 and the M3TM model 2 ^, we compare the relaxation 
rates resulting from both equations. We start with the M3TM 
equation 



dm 
dt 



R- 



Tc 



1 — m coth 



(Al) 



where we identify the Brillouin function for the case 5 = 1/2, 
B X j2 = tanh(g) with q = q l / 2 = (Tc/T e )m. Now, we use the 
identity B X j 2 ~2 /B[^ 2 smn (2<?) to write 



dm 
dt 



-R T JL 
Tc 



sinh (2q) 



>l/2 



B' 



1/2 



we multiply and divide by g/J a tp to obtain 



dm 
dt 



-R 



Tp Hat 
Tc k B T c 



2q 



sinh (2q) 



l- 



1/2, 



(A2) 



(A3) 



1 (m 2 -m 2 ) 

2 m e 



bm 



(m 2 



2m 2 



and, 



^-B\i%lm 1 (m 2 -m 2 ) 



1/2 



2 X|| 



(A5) 



(A6) 



Finally, collecting the equations JA3I > and dA6l ) altogether: 



dm ( 3R /J at \ 2T p 2q 
~dt ~ [Yk^J 37b sinh (2q) \2% 



1 m 2 

(l-—)m) (A7) 

mi 



Comparing this to the LLB equation with longitudinal re- 
laxation only and without anisotropy and external fields, we 
can write Eq. (IA7t in terms of n: 



dn X 2T e 2q 
dt ^ m e 3Tc sinh (2q) 



H e g = ya,\\H eS 



(A8) 



where H eS = & 



(1 — n 2 )n, and 




t [ps] 

Figure 6. Longitudinal relaxation calculated with M3TM and LLB 
(S=l/2) models for Nickel parameters 2 ^ and T /Tc — 0.8. 



ay = 



3R fim T p 



2y lc B T c T e 



2T e 



2q 



3Tc sinh (2q) 



(A9) 



Thus the Koopmans' M3TM equation is equivalent to the LLB 
equation with 5=1/2 and where the precessional aspects are 
not considered. The link between both of them is the identifi- 
cation 



3R flat T p 
2y^7t T e 



(A10) 



As an example we compare the result of the longitudinal re- 
laxation in a numerical experiment for both M3TM and LLB 
(5=1/2) equations. The system is put in a saturated state with 
5-/5 = 1 and we let it relax towards the equilibrium state. The 
comparison of the results for the temperature T /Tc = 0.8 are 
presented in Fig|6] 
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